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Finite Volume Procedure 
Basics 



The Finite Volume Procedure for a fluid network is an 
extension of single control volume analysis of mass and 
energy conservation in classical thermodynamics. 


Mass in 



System boundary 
Control \ \ (control surface) 


Mass out 


Control Volume Analysis 
in Classical Thermodynamics 



Finite Volume Analysis in 
Fluid Network 
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Network Definition & Main Characteristics 
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Nitrogen 70 °F 


m = 2.63 Ib/s 


* Data Structure for 
Generalized Flow 
Network 

* Pressure Based 
Finite Volume 
Method 

* Conjugate Heat 
Transfer 


Fluid Transient 



Mathematical Formulation 


Unknown Variables 

1. Pressure 

2. Flowrate 

3. Fluid Temperature 

4. Solid Temperature 

5. Mass 


Governing Equations to Solve 

1. Mass Conservation Equation 

2. Momentum Conservation Equation 

3. Energy Conservation Equation of Fluid 

4. Energy Conservation Equation of Solid 

5. Thermodynamic Equation of State 
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Program Structure 


Graphical User 


Solver & Property 
Module 


User Subroutines 



Fluid Option 

Two Thermo-dynamic Property Programs are integrated with GFSSP 

GASPAK 


GASP & WASP 


Index 

Fluid 

Index 

Fluid 

1 

Helium 

7 

Argon 

2 

Methane 

8 

Carbon dioxide 

3 

Neon 

9 

Fluorine 

4 

Nitrogen 

10 

Hydrogen 

5 

Carbon monoxide 

11 

Water 

6 

Oxygen 

12 

RP-1 


Additional Fluid can be 
added as Property Table 
Lookup 


Index 

Fluid 

Index 

Fluid 

1 

Helium 

19 

Krypton 

2 

Methane 

20 

Propane 

3 

Neon 

21 

Xenon 

4 

Nitrogen 

22 

R-11 

5 

Carbon monoxide 

23 

R-12 

6 

Oxygen 

24 

R-22 

7 

Argon 

25 

R-32 

8 

Carbon dioxide 

26 

R-123 

9 

Parahydrogen 

27 

R-124 

10 

Hydrogen 

28 

R-125 

11 

Water 

29 

R-134A 

12 

RP-1 

30 

R-152A 

13 

Isobutane 

31 

Nitrogen trifloride 

14 

Butane 

32 

Ammonia 

15 

Deuterium 

33 

Ideal gas 

16 

Ethane 

34 

Hydrogen peroxide 

17 

Ethylene 

35 

Air 

18 

Hydrogen sulfide 
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Resistance Option 



Option 

Type of Resistance 

Input Parameters 

Option 

Type of Resistance 

Input Parameters 

1 

Pipe flow 

L (in), D (in), e ID 

13 

Common fittings and valves 
(two K method) 

D(in), K v K 2 

2 

Flow-through restriction 

C L , A (in 2 ) 

14 

Pump characteristics* 

A>- % O v A (in 2 ) 

3 

Noncircular duct 

a (in), b (in) 

15 

Pump power 

P(hp), rj, A (in 2 ) 

4 

Pipe with entrance and exit 
loss 

L (in), D (in), e/D, K., K 

16 

Valve with given C, 

c,* 

5 

Thin, sharp orifice 

D \ (in), D 2 (in) 

17 

Joule-Thompson device 

L n ,V„k v ,A 

6 

Thick orifice 

L (in), D 1 (in), D 2 (in) 

18 

Control valve 

See example 12 datafile 

7 

Square reduction 

D, (in), D 2 (in) 

19 

User defined 

A (in 2 ) 

8 

Square expansion 

D, (in), D 2 (in) 

20 

Heat exchanger core 

A f (in 2 ),A (in 2 ), A (in 2 ), L 
(in), K, K 

9 

Rotating annular duct 

1- (in), r o (in), r (in), N (rpm) 

21 

Parallel tube 

L (in), D (in), f/D, n 

10 

Rotating radial duct 

L (in), D (in), N (rpm) 

22 

Compressible orifice 

C L , A (in 2 ) 

11 

Labyrinth seal 

r (in), c (in), m (in), n, a 

23 

Labyrinth seal, Egli correlation 

r (in), c (in), m (in), n, a 

12 

Flow between parallel plates 

r (in), c (in), L (in) 

24 

Fixed flow 

Flow (lb m /s), A (in 2 ) 


* Pump characteristics are expressed as Ap = Aq + B 0 iii + C 0 iir , A p - Pressure rise, lbf/ft 2 , m - Flow rate, Ibm/s . 


Additional Resistance Option can be added through User Subroutine 


Modeling Options 


Pressurization 
Heat Exchanger 
Turbopump 
Valve open/close 
Pressure Regulator 
Flow Regulator 
Cyclic Boundary 
Moving Boundary 
Fluid Mixture 
Rotation 


GRAPHICAL USER INTERFACE 



MODEL BUILDING MODEL RUNNING MODEL RESULTS 



Start I El] Introduction .ppt | g]Input_output.ppt 1 1 3 VTASC Version 2.02 - D... 6:56 PM 
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Tutorial 1 : Compressible flow in a Converging- 

Diverging Nozzle 

Detailed Schematic 

Inlet Plane Throat Exit Plane 



VTASC Model 
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GFSSP 6.05 - Tutorial 1 
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Results of Parametric Computations 






Run five cases, gradually decreasing the exit pressure. 


p . 

A exit 

(psia) 

• 

m (lbm/s) 

134 

0.279 

100 

0.329 

60 

0.336 

50 

0.337 

45 

0.337 


•How does the choked flowrate compare to the hand-calculated value of 0.327 lb m /s? 


m = A 




f 


throat inlet 


RTi 


inlet 


2 

Kr+ij 


7 + 1 


(o 


A 


= 0.19012 in 1 161.6 


Ibf 


in 2 y 


32.174 lbm (1.281) 

Ibf -s 

2.281 

^ 2 3 0.281 

1 85.83 lbi 1460°/? 

^1.281 + 1 J 


= 0.327 


Ibm 


GFSSP 6.05 - Tutorial 1 
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Tutorial 2: CHILLDOWN OF TRANSFER LINE 



Problem Considered: 

•Time dependent Pressure, Temperature and Flow Rate history during 
chilldown 


3-SHE— “ — 






Fluid Temperature 
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GFSSP 6.05 -Tutorial5 
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Tutorial -3 

Pressurization of a Propellant Tank 


Pressurant 



\ 


Propellant to Engine 


Predict the ullage conditions 
considering heat and mass 
transfer between the 
propellant and the tank wall 

Predict the propellant 
conditions leaving the tank 


Fluid: He 
P=95 psia 
7=120 °F 


GFSSP Model 
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Model Results 

Ullage and tank wall temperature history Ullage and tank bottom pressure history 



Design Project Assignment 



Air conditioning system requirements 


Requirement 

Number 

Title 

Description 

1 

Flow rate 

Each floor requires a chilled-water flow rate of 18,000 Ibm/hr. 

2 

Flow velocities 

All flow velocities are not to exceed 5 ft/s. 

3 

Floor spacing 

Each floor is 12 ft higher than the one below. 

4 

Supply and return pipe sizing 

The supply and return pipes must all be the same diameter (Schedule 80). 
The total length of supply pipe running from the exit of the bottom floor, 

5 

Supply pipe length 

through the pump and chiller, and up to the entrance of the top floor is 300 ft 
(See Fig. 11). 

6 

Floor pipe sizing 

The floor pipes must all be the same diameter (Schedule 80). 

7 

Floor pipe length 

Each floor must have 100 ft of pipe. 

8 

Floor valves 

Each floor must have two gate valves. 

9 

AHU pressure drop 

The pressure drop across each AHU must be 50 psi. 

10 

Chiller pressure drop 

The pressure drop across the chiller must be 100 psi. 
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GFSSP Model 
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Project design results 



GFSSP Value 

Requirements Value 

Notes 

Supply and return pipe size 
(Schedule 80) 

3.826 in 



Floor pipe size 
(Schedule 80) 

1.939 in 



Pressure rise across pump 

154.0 psi 



Pump power required 

23.5 hp 
(17.5 kW) 


At 55% efficiency 

Cost of operation (per 24 hr) 

$42 


$0.10 per kW-hrfor 24 hr 

Mass flow rate in pipe 1011 
(4th floor) 

18,036 Ibm/hr 

18,000 Ibm/hr 

% Difference: 0.2% 

Mass flow rate in pipe 3031 
(3rd floor) 

18,018 Ibm/hr 

18,000 Ibm/hr 

% Difference: 0.1% 

Mass flow rate in pipe 3738 
(2nd floor) 

18,018 Ibm/hr 

18,000 Ibm/hr 

% Difference: 0.1% 

Mass flow rate in pipe 4344 
(1st floor) 

18,040 Ibm/hr 

18,000 Ibm/hr 

% Difference: 0.22% 

Mass flow rate in 300-ft 
supply pipe 

72,108 Ibm/hr 

72,000 Ibm/hr 

% Difference: 0.15% 

AHU 4 pressure drop 

49.96 psi 

50 psi 

% Difference: 0.08% 

AHU 3 pressure drop 

49.86 psi 

50 psi 

% Difference: 0.28% 

AHU 2 pressure drop 

49.86 psi 

50 psi 

% Difference: 0.28% 

AHU 1 pressure drop 

49.96 psi 

50 psi 

% Difference: 0.08% 

Chiller pressure drop 

100.1 psi 

100 psi 

% Difference: 0.1% 

Fluid velocity in pipe 1011 
(4th floor) 

3.917 ft/s 

<5 ft/s 

Meets requirement 

Fluid velocity in pipe 3031 
(3rd floor) 

3.913 ft/s 

<5 ft/s 

Meets requirement 

Fluid velocity in pipe 3738 
(2nd floor) 

3.913 ft/s 

<5 ft/s 

Meets requirement 

Fluid velocity in pipe 4344 
(1st floor) 

3.917 ft/s 

<5 ft/s 

Meets requirement 

Velocity in 300-ft supply line 

4.022 ft/s 

<5 ft/s 

Meets requirement 


Summary 






• This paper demonstrates how a system level, user-friendly 
network flow analysis code can be integrated in a senior level 
thermal design class. 

• The intent was to introduce a state-of-the-art computational 
tool to perform a real world technical task 

• The introduction of GFSSP was done through lectures, 
tutorials, and a senior design project. 

• The authors' experience of using GFSSP in the class was very 
positive from students'feedback 

• GFSSP is available free of cost to all universities in the United 
States from MSFC's Technology Transfer Office 
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